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The paper by Mayoh and Garcia-Garcia arXiv:1412.0029irl] is entitled ’’Can disorder re¬ 
ally enhance superconductivity?”. In our opinion, the answer given by the authors is not 
satisfactory. We present the alternative picture. 


The paper by Mayoh and Garcia-Garcia [I] is en¬ 
titled ’’Can disorder really enhance superconductiv¬ 
ity?” . In our opinion, the answer given by the au¬ 
thors is not satisfactory and we present the alterna¬ 
tive picture. 

Mean field solution. A basis for description of the 
spatially inhomogeneous superconductivity is given 
by the Gor’kov equation for the order parameter 
A(r) 

A(r) = J K{r,r')A{r')d‘^r' (1) 

with the kernel K{r,r') satisfying the sum rule [2] 

J K{r,r')d‘^r' = (r)ln , (2) 

where g is the Cooper interaction constant, is 

the local density of states at the Fermi level, ojq is a 
cut-off frequency, and d is a dimension of space. 

The Anderson theorem [3] follows from Eq.l un¬ 
der assumption of a self-averaging order parameter, 
when A(r) and K(r,r') can be independently aver¬ 
aged over disorder. Since (A(r)) does not depend on 
r due to the spatial uniformity in average, the use of 
the sum rule (2) gives 

{A)=g{uF)\nhlp^{A), (3) 


one corresponds to the moderate variation of A(r) 
in space, so that it remains more or less of the same 
order in the whole space. The corresponding Tc is 
given by the formula 


STc _ 1 f ^ i/oUi{r) + iyi{r)'^ 

T,o " ALd i ^ ^ 


(4) 


where i^i(r) is a deviation of the local density of 
states vpir) from its unperturbed value vq, \ = gi/Q 
is the dimensionless coupling constant, is the 
transition temperature in the absence of disorder, 
is a volume for one impurity, and integration is 
carried out over a vicinity of the single point de¬ 
fect. The linear in viir) term exactly corresponds 
to the Anderson theorem and relates the change in 
Tc with the change of the average density of states. 
Generally, vi{r) is comparable with vq and already 
Eq.4 predicts a possibility of essential violation of 
the Anderson theorem. It is related with the fact 
that the initially uniform order parameter is influ¬ 
enced by point defects and can increase or decrease 
in their vicinity. 

More essential deviations from the Anderson the¬ 
orem arouse in the second regime, when the order 
parameter is mainly localized at the small number 
of ’’resonant” impurities producing the quasi-local 
states near the Eermi level. The corresponding esti¬ 
mate for Tc [4] 


and the critical temperature Tc is given by the BCS 
formula, which contains the average density of states 
{vp)- If the latter is maintained fixed, Tc is not 
changed by disorder. However, self-averaging does 
not hold in the general case, and deviations from 
the Anderson theorem arise. 

Equation (1) can be accurately solved for a small 
concentration of the point-like impurities [4]. This 
solution shows possibility of two regimes. The first 


Tc - ga-'^ - ATo (5) 

(a is the lattice spacing, Eq is a scale of the order 
of the Eermi energy or the bandwidth) is valid for 
small g and saturates by a quantity of the order ujq, 
when g increases. 

Results (4) and (5), obtained for a small cocen- 
tration c of impurities, can be qualitatively extrap¬ 
olated into the c ~ 1 region. It allows to give the 
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adequate answer to the question in the title of the 
paper: 

(a) Disorder can enhance Tc by a trivial reason, 
due to increase of the average density of states. 

(b) If the average density of states remains fixed, 
then Tc can be changed due to deviations from the 
Anderson theorem. This effect is always positive in 
the framework of formula (4). 

(c) There is a possibility of the catastrophic in¬ 
crease of Tc due to resonances on the quasi-local 
levels, though this regime is affected by fluctuations 
(see below). 


Role of multifractality. Recently there have been 
claims [SI 1 [7] that a great increase of Tc is pos¬ 
sible in the vicinity of the Anderson transition due 
to multifractality of wave functions. The analysis of 
[SJE] is based on equation 


A(e) = - 


A f /(e,e')A(c') 


x/e'" + A(eO 

— UJ{} 

with the kernel of the form 

7(e,e') = 


tanh 


Ve'" + A(e')' 


2T 


En 


j de' 
( 6 ) 

(7) 


motivated by multifractal properties of wave func¬ 
tions. Tc is determined by equation ( 6 ) linearized in 
A: accepting A to be a function e/T and dimension- 
alizing the integral, one has for infinite ujq |Sli 

Tc - . ( 8 ) 


The singular limit of small 7 was analyzed in [I] and 
lead to result 


Tc ^ Eq 




(9) 


which reproduces ( 8 ) in the limit ujq 00 , the BCS 
result in the limit 7—^0 and describing saturation 
of ( 8 ) by a value ^ wq for large A. The claim of [1] 
that ujQ is always less than Eq is incorrect: Eq can 
be small in semimetals and narrow band materials. 

The linearized form of equation ( 6 ) is not equiva¬ 
lent to ( 1 ) (see discussion in [ 6 ]) and based on two as¬ 
sumptions: (i) truncation of the Hamiltonian in the 
BCS spirit, and (ii) averaging of the kernel indepen¬ 
dently of aB Both approximations are uncontrol- 

^ It means a self-averaging assumption, but in the modified 
form: instead of the usual equation (A) = (K){A) one uses 
(SA) = {SKS~^)(SA) where 5 is a certain operator. It can 
be interpreted in the variational spirit, considering 5 as a kind 
of the trial function. 


lable. The partial justification of (i) was suggested 
in [I]: truncation of the Hamiltonian is rigorous for 
a pure superconductor and should be valid approx¬ 
imately in the case of weak spatial inhomogeneity. 
One can partially excuse assumption (ii), using some 
effective exponent in the capacity of 7 . Indeed, the 
kernel I(e,e') is determined by a matrix element 
J c?^r|' 0 (e, r)p|' 0 (e',r)p whose estimation gives dif¬ 
ferent values of 7 for the ’’usual” and ’’typical” aver¬ 
aging 0; this uncertainty is aggravated, if averaging 
is made with the weight A(r). 

The first argument explains why the strong local¬ 
ization regime (corresponding to (5)) does not con¬ 
tain in ( 6 ) for I(e, e') = const. The second argu¬ 
ment shows that a difference between (5) and ( 8 ) 
arises on the level, which is not controllable in the 
approach based on Eq. 6 . In fact, the difference be¬ 
tween (5) and ( 8 ) is not quantitative, but qualitative. 
Result (5) is not restricted by a vicinity of the mo¬ 
bility edge and remains approximately the same for 
any energy inside the band; correspondingly, it has 
no relation to multifractality. Nevertheless, the real 
physical mechanism is the same for results (5) and 
( 8 ) and related with resonances at quasi-local levels 
[3]. Indeed, if the local density of states v{e, r) is con¬ 
sidered as a smooth function of e, then its variation 
is finite: it corresponds to shifts of the whole band 
by a value W or —W for the Anderson model with 
distribution of site energies in the interval (—IT, IT). 
Unbounded fluctuations of z^(e, r), arising in the con¬ 
text of multifractality, are necessarily related with a 
partial discretization of the spectrum due to a pres¬ 
ence of quasi-local levels. If the usual value of 7 is 
exploited in (7), then the Tc value given by (5) is 
greater than ( 8 ); it means that the Cooper instabil¬ 
ity occurs at configurations, which are governed by 
individual peaks and not fractal clasters | 3 ]. 

Weak multifractality considered in [I] is practi¬ 
cally actual only for the 2D case in the regime of 
weak disorder. However, this regime is described 
by formula (4), which can be obtained in this case 
making two iterations of the Gor’kov equation and 
exploiting the sum rule (2). In this case one can 
integrate along all the system and use its whole vol¬ 
ume in the capacity of L‘^: dividing of disorder into 
separate ’’impurities” is not necessary. There is no 
need to use the approximate equation ( 6 ), when the 
accurate result is available. By the way, in the frame¬ 
work of (4) the order parameter is proportional to 
i'p{r) and the logarithmically normal distribution 
for A(r) [T] follows trivially in the weak multifrac- 
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tality regime. 

It is clear from this consideration, that multifrac- 
tality is not a direct cause for the increase of Tc and 
results of kind (8) are related with the more univer¬ 
sal mechanism. 

Role of fluctuations. Equation (5) is a result of 
the mean field theory. The corresponding configura¬ 
tion of the order parameter is a uniform background 
Aq with abrupt peaks at resonant impurities, whose 
concentration is of the order Tc/Eq. The order pa¬ 
rameter can be considered as positiv^, and so its 
phase is the same in the whole volume. When we 
come to a fluctuational description, the modulus of 
the order parameter remains practically unchanged, 
while the essential phase fluctuations arise. If the 
uniform background is neglected, then the system is 
divided into practically independent superconduct¬ 
ing ” drops”, whose phases are fluctuating freely and 
destroy the macroscopical coherence of the super¬ 
conducting state. If the uniform contribution Aq 
is taken into account, the Josephson coupling be¬ 
tween drops arises and their phases become corre¬ 
lated. The accurate fluctuational analysis of such 
a system is nontrivial, but the general character of 
results is the same as for the granular superconduc¬ 
tors [8]. If the ratio Tc/E q is not too small, then 
the resonant impurities are close to each other and 
their Josephson interaction is strong enough for sta¬ 
bilization of the mean-field solution at practically 
the same Tc value. Contrary, if TdE q is sufficiently 
small, then the Josephson coupling between drops is 
weak and fluctuations destroy superconductivity at 
temperatures close to the mean-field Tc value. How¬ 
ever, decreasing of temperature stimulates the grow¬ 
ing of tails of the localized solutions [5]; the Joseph¬ 
son coupling between drops increases and stabilizes 
superconductivity before Tcq is reached. Hence, fluc¬ 
tuations suppress Tc in comparison with its mean- 
field value, but do not eliminate enhancement of Tc 
completely. 

Analogous arguments can be given for configura¬ 
tions corresponding to (8), where a fraction of the 
superconducting phase is estimated as {Tc/Eod [6]. 
However, such configurations are not actual, since 
the Cooper instablity corresponds to (5). 

The paper [1] suggests another way to deal with 
fluctuations. Firstly, solution of (6) for T = 0 is 


^ In the absence of magnetic effects, the kernel K{r,r') is 
positive, and the Cooper instability corresponds to the node¬ 
less eigenfunction. 


found, giving the spatially inhomogeneous order pa¬ 
rameter Ao(r). Secondly, the field Tc{r) of ’’local 
Tc” values is introduced, such as Tc{r) oc Ao(r).0 
Finally, the global Tc is defined as a percolation 
threshold in the field Tc{r). Such percolation pic¬ 
ture has a sense for certain conditions m. but it is 
not the case for weak spatial inhomogeneity. 

Indeed, the Gor’kov equation (I) defines the spa¬ 
tially inhomogeneous configuration A('r), which ap¬ 
pears at a certain temperature. This temperature, 
by definition, is a final result for Tc'. there is no need 
to consider ’’local Tc” or ’’percolation”. Of course, 
one should attend for insignificance of fluctuations, 
but this condition is rather weak. The Ginzburg 
number is incredibly small for a pure superconduc¬ 
tor; it increases gradually with increasing of spatial 
inhomogeneity, but it is possible to reach rather large 
values of ratio lS.max /^min before this number will 
approach unity and fluctuations will become essen¬ 
tial. Surely, no percolation is necessary for weak 
spatial inhomogeneity. The percolation picture be¬ 
comes reasonable for the Ginzburg number of the 
order of unity, when the mean-held estimate of Tc is 
poorly dehned and the use of percolation allows to 
rehne it. 

Role of interaction. The BGS constant A corre¬ 
sponds to some effective interaction. In a more de¬ 
tailed description it is combined from the electron- 
phonon coupling and the Goulomb pseudopotential. 
The latter is known to increase in the presence of 
disorder and it is the main cause for Tc degradation 
[min]. This effect was not considered explicitly in 
[DllElle], but is surely essential in discussion of the 
experimental situation. It is the main reason why 
enhancement of Tc is a rare thing in reality. 

In paper [7] the result analogous to (8) is obtained 
in the framework of the Finkelstein renormalization 
group approach. However, it is completely different 
from [gig in the initial assumptions and the dis¬ 
cussed physical mechanism. The Gooper constant g 
is kept hxed in [Si, but an attempt is made to ad¬ 
vance beyond applicability of the Anderson theorem. 
Gontrary, the authors of [7] consider renormalization 
of g by disorder, while self-averaging is accepted for 
granted. By the latter reason, the strongly localized 
regime was not accessible in this approach, while the 
obtained effect is lesser than (5). The questions also 
arise, how results of agree with the usual picture 

® In fact, this relation is violated due to the presence of 
scale Tco- 
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of the Coulomb pseudopotential enhancement. 

In conclusion, disorder can enhance T^. (a) due 
to increase of the average density of states; (b) due 
to deviations from the Anderson theorem; (c) due 
to resonances on the quasi-local levels. The latter 
regime is affected by fluctuations, and in some cases 
can be essentially suppressed. Practically, enhance¬ 
ment of Tc by disorder is a rare thing due to the 
increase of the Coulomb pseudopotential. 

The author is indebted to M.V.Sadovskii for the 
discussion of paper [1] . 
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Reply to arXiv: 1502.06282 

As a reaction to the preceding overview m, Mayoh 
and Garcia-Garcia have submitted the comment m- 
This comment has a form of the personal attack and 
contains the whole series of untrue statements. We 
give brief remarks, in order to reveal these statements. 

1. Our mean field results for a superconductor 
with a small concentration of the point-like impu¬ 
rities [4] are based on the accurate solution of the 
Gor’kov equation. These results still persists, even 
if somebody does not like them. Gontrary, the ap¬ 
proach of [T] involves uncontrollable approximations 

m- 

2. It is repeatedly stated in [14] that we pre¬ 
dict Tc ^ bOOOK for the superconducting transition 
temperature, which ’’has never, and very likely will 
never, been found numerically or experimentally”. 
In fact, it was clearly indicated (see the text after 
Eq.2 in |4]) that Tc is bounded by the quantity ujq/tt 
{ ojQ is the cut-off frequency), which for the phonon 
mechanism corresponds to values already observed 
for oxide superconductors {Tc ^ 150K, ujq ^ 400K). 
A possibility to describe the latter in terms of our 
model was discussed in Sec.7 of |4]. 

3. The authors of m write on ’’unjustified use 
of the mean field approximation”. In fact, insuffi¬ 
ciency of the mean field approach was indicated in 

and the role of fluctuations was extensively 
discussed. We cannot help, if our arguments are ig¬ 
nored. 

4. According to [14], our formula (4) in [13] is re¬ 

lated with ’’poorly defined” quantities and unclear 
conditions of applicability. In fact, all explanations 
were given in mm and we can repeat those of them 
that are questioned in [14]: is a volume for one 

impurity, i.e. the whole volume divided for a num¬ 
ber of impurities; the quantity i^i{r) is the difference 
iyp{r) — i/Q, where vpir) is the local density of states, 
defined in a standard manner (see Eq.8 in [4]), and 
z/Q is its unperturbed value. For a small amplitude 
of disorder (|z/i(r)| <C zzq) the indicated formula can 
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be obtained by two iterations of the Gor’kov equa¬ 
tion m- It remains valid for a small concentration 
of strong impurities, when |i^i(r)| ^ i^q [1]. It is 
not restricted by any assumption on the correlation 
between impurities. 

5. According to [14], ’’Suslov claims ... that we 
do not discuss phase-fluctuations”. There is no such 
statement in [T3]; our comments on the role of inter¬ 
action are brief and contain no criticism of [T]. 

6 . According to [14j . there is a wrong statement in 
m- ” Suslov also states that we claim Eq is always 
larger than e£)”. One can compare it with the text in 
[Tj: ”We believe that this is necessary since eo ^ Eq 
so it is inconsistent to take the Debye energy to in¬ 
finity while keeping Eq finite”. This argument was 
used in [T] to reject the result Tc ^ EqX^^'^ appear¬ 
ing in preceding publications: ’’This is an expression 
that, we are at pains to stress, is not recovered in 
our formalism” M- In fact, this result follows from 
Eq.22 of [1] in the limit eo —t oo, independently of 
the desire of authors. 

7. According to [M], ’’Suslov claims that our re¬ 
sults are only valid in small and strictly two dimen¬ 
sional systems”, while they are valid also for thin 
films and two-dimensional systems with spin-orbit 
interactions. In fact, it is written in m that ’’weak 
multifractality considered in [1] is practically actual 
only for the 2D case in the regime of weak disorder”, 
and there is no rejection of the indicated additional 
applications. 

8 . Our main criticism of [T] refers to the use of the 
percolation picture uni beyond the limits of its appli¬ 
cability. Indeed, the case of the ’’moderate spatial 
inhomogeneity” is described by the Gor’kov equa¬ 
tion, which has solution A(r) arising at a certain 
temperature. This temperature, by definition, is a 
final result for Tc'. there is no need to consider ’’lo¬ 
cal Tc” or ’’percolation”. Such situation persists, till 
fluctuations are insignificant and the mean field esti¬ 
mation of Tc is well-defined. The percolation picture 
becomes reasonable for ’’very strong spatial inhomo¬ 
geneity” , when the Ginzburg number becomes of the 
order of unity; in this case the mean-field estimate of 
Tc is poorly defined and the use of percolation allows 
to improve it. 

The authors of [T] try to disprove our arguments 
basing on the difference between ’’weak inhomogene¬ 
ity” and ’’weak multifractality”. This attempt is not 
successful due to the facts: 


(a) The suggested in [T] partial justification of 
Eq.6 in m refers to ” weak inhomogeneity” and not 
to ’’weak multifractality”. 

(b) It is clear from [T3] that for applicability of 
percolation, the ratio A^axl^mm should be greater 
than some large parameter; here A^ax and Amin are 
typical (not exclusive) values. The distribution for 
A(r’) is presented in Fig.4 of [T]. This distribution 
can be cut-off on both sides, since its tails corre¬ 
spond to local perturbations, which have no conse¬ 
quences for global superconductivity. After that the 
ratio Amax/Amin is typically of the order of unity: 
it corresponds to ” moderate inhomogeneity” and ap¬ 
plicability of the Gor’kov equation. 

(c) The situation can be discussed constructively. 
Practically, ’’weak multifractality” corresponds to 
the weakly disordered 2D case, which is described 
by formula (4) in |13j . The corresponding order pa¬ 
rameter is either slightly perturbed (for weak im¬ 
purities), or its perturbations are local (for a small 
concentration of strong impurities). In both cases, 
there are no problems with fluctuations, and hence 
there is no place for percolation. 

The authors of m write correctly that ’’accord¬ 
ing to Suslov our percolation analysis is intended to 
describe fluctuations”: it is right, in spite of their 
objections. They are wrong in ascribing to us an 
idea that ’’the value of Tc resulting from percolation 
is similar to that obtained by averaging over Tcir)". 

9. It is stated in m that ”in weakly coupled su¬ 
perconductor A <C 1 so Tc/To is always small and the 
global critical temperature must be necessary zero”. 
This statement reveals a complete misunderstanding 
of our arguments. We say that, for not very small ra¬ 
tio Tc/Eq, the average distance a (Eq/Tc)^^^ between 
resonant impurities becomes comparable with a and 
the background value Aq of the order parameter be¬ 
comes comparable with its resonant peaks. Then our 
’’strongly localized regime” becomes ’’moderately lo¬ 
calized” and has no problems with fluctuations. In 
any case, we do not restrict ourselves by weak cou¬ 
pling, and in the worst situation Tc falls to Tco and 
not to zero. 
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